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Quantum spin chains with site dissipation 
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We use Monte Carlo simulations to study chains of Ising- and XY-spins with dissipation 
coupling to the site variables. The phase diagram and critical exponents of the dissipative Ising 
chain in a transverse magnetic field have been computed previously, and here we consider a 
universal ratio of susceptibilities. We furthermore present the phase diagram and exponents 
of the dissipative XY-chain, which exhibits a second order phase transition. All our results 
compare well with the predictions from a dissipative field theory. 
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Dissipation and decoherence in mesoscopic quantum 
systems have important practical implications, from sta- 
bilizing superconductivity in granular materials, to the 
loss of information stored in qubits. The macroscopic de- 
gree of freedom can often be described by a spin variable 
of Ising or XY symmetry. They might represent the low- 
est energy states in a double well potential, or the phase 
of the superconducting order parameter on some grain 
or island. Following Caldeira and Leggett"'^ one can in- 
troduce dissipation in such a system by coupling this 
spin variable linearly in the coordinates to an infinite set 
of harmonic oscillators which represent the environment. 
The spectral distribution of these harmonic oscillators 
defines the type of dissipation and we will consider a so- 
called Ohmic heat bath which produces linear friction. 
Upon integrating out the heat bath degrees of freedom 
one obtains an imaginary time effective action with long 
ranged interactions in time. 

The simplest of these models, the dissipative two-state 
(Ising) system, has been studied extensively, both the- 
oretically and by numerical simulations.^^'* It exhibits 
a dynamical phase transition to a localized state at a 
critical value of the dissipation strength. We have re- 
cently studied what happens when infinitely many such 
two-states systems are spatially coupled.^ We observed 
a quantum phase transition in a universality class differ- 
ent from either the single site case or the non-dissipative 
system. 

The second class of models we shall study involve 'spin' 
degrees of freedom with a global U(l) rotation (XY) sym- 
metry. A single XY-spin coupled to dissipation describes 
the physics of the so-called single-electron box.^'^ The 
integer variable conjugate to the compact phase repre- 
sented by the spin corresponds to the number of excess 
charges on the small metallic island, which is connected 
to an outside lead by a tunnel junction. XY spins on 
higher dimensional lattices describe arrays of Joseph- 
son junctions. In resistively shunted junctions, the dis- 
sipation couples to the phase difference and the phases 
are usually treated as non-compact variables. Such sys- 
tems exhibit a dissipation driven superconductor-to- 
metal transition'^^^* because the dissipation suppresses 



phase slips and can restore superconductivity. Very re- 
cently, an XY-chain with site dissipation was studied in 
connection with the superconductor-to-metal transition 
in nanowires.*^ 

In this work, we shall consider chains of Ising- and XY- 
spins with Ohmic dissipation coupling to the site vari- 
ables. The discretized action for a chain of n-component 
spins with on-site dissipation reads 
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, n and cr^ = 1 on every site. An 
Ising chain (in a transverse magnetic field) corresponds 
to n=l, a chain of XY-spins to n = 2, and so on. The 
number of lattice sites in the spatial and imaginary time 
directions are Nx and Nt, respectively, and we apply pe- 
riodic boundary conditions in both directions. The long- 
ranged term containing the dimensionless variable a de- 
scribes the Ohmic dissipation coupling to the site vari- 
able CT. 

The nature of the phase transition and its universal 
properties have recently been discussed by Pankov et 
al.,^^ who anayzed the dissipative n-component 0"* field 
theory (the chain (1) corresponds to spatial dimension 
d=l) 
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They predict a second order phase transition for a > 0, 
whose properties should not depend on the dissipation 
strength. The upper critical dimension of the field theory 
(2) is 2, independent of n, and an expansion to second 
order in e = 2 — d yields the following expressions for the 



2 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



critical exponents^^' 
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An additional theoretical prediction concerns the 
Fourier transform of the spin-spin correlation function, 
or susceptibility x(A;,ia;„). Here, k denotes a wave vec- 
tor and Un = 2-Kn/NT is the Matsubara frequency (the 
absolute temperature is T = h/{kBN-r)). The strong hy- 
perscaling properties of the field theory in Eq. (2) have 
interesting consequences for its universal properties at 
the quantum critical point. In particular, because the 
last |a;| term in Eq. (2) does not renormalize, there is no 
non-universal cut-off dependence to the overall scale of 
X{k,uin), which thence obeys the scaling form 
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where only the number ci is non-universal, and all other 
aspects of the scaling function $ are independent of cut- 
off scale parameters. In the present situation, we are not 
simulating the field theory in Eq. (2), but rather the 
lattice Ising model in Eq. (1), and so we need a relation 
between the overall scale of 4) and the Ising spins. We 
sidestep the ignorance of the latter scale by computing 
the ratio 



T(n) = 



(5) 



x(o,o) ■ 

In the limit of an infinite spatial length (1 ^ A^.^ ^ ^x)i 
Eq. (4) implies the remarkable result that the function 
T(n) is a completely universal function of the integer n. 
We can also determine the limiting behavior of this func- 
tion for small and large n. At low frequencies, fujj <c fcsT, 
we expect dissipative behavior with a linear \ujn\ depen- 
dence of the dynamic susceptibility. At large JvjJ ^ AsT, 
scaling implies^^ x(0, ~ l/|a;„|(2-»))A 
Hence we have 



for 1 <C n < A^T- 



T(n) = i / n r 

^ ' 1 1 - n/B for n < 1, 



(6) 



where A and B are universal numbers. 

We will first test the universality of the function T for 
Ising spins and then present the phase diagram and criti- 
cal exponents for the XY-chain. In the Monte Carlo sim- 
ulations, we use a variant^^ of the Swendsen-Wang algo- 
rithm^"'' which builds the clusters in a time 0{Nt log Nt) 
despite long ranged interactions. 

Since n <C 1 is not accessible in imaginary time, we 
consider the case 1 <C <C Nr and compute the am- 
plitude A for weak, intermediate and strong dissipation, 
setting Kr = — ^ ln[tanh(l)] = 0.136 as in our recent 
study of the Ising chain. ^ For this value of Kr, the phase 
transition for a single dissipative Ising spin occurs at 
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Fig. 1. T(n) for a = and a = 0.3. The different curves corre- 
spond to different Nt in the range 1 <C A'^t ^ {a = 0) and 
1 <C Nt <ti N^ (a = 0.3) respectively. The data are plotted in 
the range 1 < n < Nt/2. 
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Fig. 2. T(n) for a = 0.3 and ^ ^ ^ 

sets corresponding to different A'^r are plotted in different colors 
for better clarity. The dotted line shows the best fit to A/n. 



oic = 0.625. The universal function T, however, should 
not depend on the choice of Kt- 

For a (non-dissipative 2D Ising model), the am- 
plitude A can be calculated exactly. One finds for large 
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which we used to test the numerical simulations. (Note 
that the large n behavior ^ n^^/* applies only for a = 0; 
any non-zero a is expected to crossover into the behavior 
in Eq. (6).) In this case the dynamical critical exponent is 
2 = 1, so we plot in Fig. 1 the function T(n), computed 
at the critical inter-site coupling for different lat- 
tices with 1 <^ Nr <^ Nx- The behavior predicted from 
the Onsager solution is nicely reproduced in the scaling 
limit 1 <C n ^ Nr ■ Also shown in Fig. 1 are the data for 
a = 0.3. Since z w 2 for a > 0, we chose lattices with 
1 <^ Nt ^ N^. A power-law decay proportional to 1/n 
is clearly visible, although the curves bend away slightly 
for n outsite the scaling region 1 <^ n <^ Nr- For the pur- 
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Fig. 3. Data collapse of Binder cumulants at the critical point for 
a = 0.3 and z = 1.97. 




Fig. 4. Phase diagram of the dissipative XY-spin chain in the 
space of nearest neighbor coupling and dissipation strength 
a {Kr = 0.1). 



pose of extracting the amplitude A, we define the scaling 
region as <n< j^-^- Fig. 2 shows the data for 

a = 0.3 restricted to these values of n. A fit to the pre- 
dicted power-law (6) yields A = 0.0729. Repeating the 
same analysis for a = 0.1. wc get A = 0.0726. For very 
strong dissipation (a = 0.5) we see stronger corrections 
to scaling and n ~ lies clearly outside the scaling 

region. Larger lattices would be needed to extract an ac- 
curate value of A. Nevertheless, even our data at a = 0.5 
are roughly consistent with the asymptotic form 

0.073 



T(n) 



(8) 



Having confirmed the predictions for T, we use quan- 
tum Monte Carlo simulations to determine the phase di- 
agram and critical exponents of the dissipative XY-spin 
chain following the analogous investigation of the Ising- 
spin case in Ref. 5. Assuming a second order phase tran- 
sition, wc use a scaling procedure which was proposed 
for the study of spin glasses. For a given value of a 
{Kr =0.1 is kept fixed in this study) we determine the 
critical spatial coupling and the dynamical critical 
exponent z self-consistently from the data collapse of the 
Binder cumulant ratio B = 1 — (m*)/(3(m^)^), where m 
is the magnetization. Near the critical point, B scales as 

BiN.M=^s{^,^). (9) 

At the critical point, the correlation length ^ diverges, 
and the Binder cumulants collapse onto a universal func- 
tion of N-r / {NxY ■ The excellent data collapse shown in 
Fig. 3 validates the scaling Ansatz in Eq. (9) and con- 
firms the second order transition which was predicted in 
Pankov eA, al}^ The data collapse yields = 0.92132(2) 
and z = 1.97(3). 

Repeating this analysis for several valuc;s of a, we map 
out the phase diagram shown in Fig. 4. In contrast to 
the Ising case^"^ there exists no ordered state at = 0. 
As mentioned in the introduction, the single site case 
corresponds to the imaginary time effective action of a 
single-electron box,^ which was recently studied using 
extensive Monte Carlo simulations.^ 

Next, we compute the critical exponents r} and v. The 



order parameter correlation function C{x, t) in an infi- 
nite system and at the critical coupling if^ scales as 

C(.x,t) - x-^'+'^-^^-gir/x") = t-'^'-+'^-^^'' g{x' /t). (10) 

The equal-time correlation function at the critical point 
should therefore decay asymptotically as C{x,Q) ~ 
xhe corrections to scaling are larger than in 
the Ising case^ and we cannot directly fit a periodic ver- 
sion of this power-law to the correlation function. In- 
stead we consider the correlations between the most dis- 
tant sites for different system sizes. In Fig. 5 we plot 
{(Ja{Q,T)(Ja{Nx/2,T)) as a function of Nx- A power-law 
fit to these data yields z + r] = 1.985(20). 
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Fig. 5. Equal-time correlation ((To(0, r)<Ta(^x/2, r)} at the criti- 
cal point for a = 0.3. The time direction of the lattice is scaled 
as Nt- = lON^. The curve shows a power- law fit to the data. 



In order to determine the critical exponent v, we set 
Nr — 5{Nxy to reduce the two-parameter scaling to a 
one-parameter one. This allows us to employ the proce- 
dure outlined in Refs. 18 and 5. Using the definition 

[m"] = (m"S,)K= - {mnK^J^.)K^^, (11) 
where m denotes the magnetization and E^, the spatial 
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Fig. 6. Plot of 21n[m'^] — ln[m*] ([m"] is defined in (11)) as a 
function of InA'a; at the critical point (a = 0.3, Kx = 0.1). We 
scaled the imaginary time dimension as Nt = 5{Nx)'' with z=2 
and z=1.97 (not shown). The slope of the fitted line gives l/i/. 



Table I. Comparison of exponents obtained by our simulations 
and the values predicted by the e-expansion in order O(e^). The 
values for the Ising chain are from Ref. 5. 
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0.638(3) 


0.633 


0.689(6) 


0.663 
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0.015(20) 
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0.021 
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1.985(15) 


1.980 
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1.979 
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2.00(1) 
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1.985(20) 
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coupling energy, one finds from finite size scaling consid- 
erations 

2 ln[m2] - ln[m''] = const +^\nN^. (12) 

The left hand side of Eq. (12) is plotted in Fig. 6 for 
a = 0.3 and Kt = 0.1. From the slope of a linear fit to 
these data one obtains v. The major source of error is 
the uncertainty on 2; = 1.97(3). Hence, we determinded 



the slopes for z = 1.97 and z = 2 and find v = 0.689(6). 

Having obtained the critical exponents for a = 0.3, we 
would like to compare our result with the predicitions 
from the dissipative n-component 0"* field theory (2) . One 
dimensional chains correspond to e = 1. In Table I we list 
the numerically determined exponents for Ising {n — 1) 
and XY-chains {n — 2) and the values from the two-loop 
e-expansion (3). The agreement is remarkably good. In 
the Ising case we found that the exponents do not depend 
on the dissipation strength and we expect the same for 
the XY-chain as well. 
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